An accurate formula for the period of a simple pendulum 
oscillating beyond the small angle regime 

F. M. S. LimJ] 

Instituto de Fisica, Universidade de Brasilia, 
P. 0. Box 04455, 70919-970, Brasilia-DF, Brazil 

P. Arun 

Department of Physics and Electronics, SGTB Khalsa College, 
University of Nova Delhi, Delhi 110 007, India 

Abstract 

A simple approximate formula is derived for the dependence of the period of a simple pendulum 
on the amplitude. The approximate is more accurate than other simple formulas. Good agreement 
with experimental data is verified. 
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I. INTRODUCTION 



The periodic motion exhibited by a simple pendulum is harmonic only for small angle 
oscillations.— Beyond this limit, the equation of motion is nonlinear. Although an integral 
formula exists for the period of the nonlinear pendulum, it is usually not discussed in in- 
troductory physics classes because it is not possible to evaluate the integral exactly.^ For 
this reason almost all introductory physics textbooks and lab manuals discuss only small 
angle oscillations for which the approximation sin ^ ^ is valid. The linearized equation has 
a simple exact solution, whose derivation can be easily understood by first-year students.— 
This linearization has bothered us since our own undergraduate days because the amplitude 
needs to be less than 7° if an error less than 0.1% (the typical experimental error obtained 
with a stopwatch) is desired. Measurements in undergraduate labs rarely have such small 
amplitudes,— and interested students sometimes ask for a relation that can describe the 
increase of the period observed for large amplitudes.— 

The restriction to small angle oscillations hinders the understanding of real-world be- 
havior because the pendulum isochronism observed in the small angle regime vanishes for 
increasing amplitudes. This restriction is also unnecessary because millisecond precision in 
measurements of the period is easily obtained with current technology .^1^1^1^ For instance, 
an experimental error of the order of 0.1% or less is typically obtained with a one meter- 
long pendulum, and thus accurate experimental studies of the dependence of the period on 
amplitude are possible even in introductory physics labs.-^ 

In this paper we derive a simple and accurate formula for the period of a pendulum 
oscillating beyond the small angle regime. The deviation from the exact results are of the 
same order of the experimental error. 

II. APPROXIMATION 

An ideal simple pendulum consists of a particle of mass m suspended by a massless rigid 
rod of length L that is fixed at the upper end such that the particle moves in a vertical 
circle. This simple mechanical system oscillates with a symmetric restoring force (in the 
absence of dissipative forces) due to gravity, as illustrated in Fig. 1. Its equation of motion 
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is given bj*^ 

-^ + fsm» = 0. (1) 

where 6 is the angular displacement in radians = at the equilibrium position) and g 
is the local acceleration of gravity. For any given initial condition, the exact solution can 
only be obtained numerically (with arbitrary accuracy). For small angle oscillations, the 
approximation sin 6' ^ 6' is valid and Eq. (1) becomes a linear differential equation analogous 
to the one for the simple harmonic oscillator. In this regime the pendulum oscillates with 



a period Tq = 2tx^ Lj a well-known textbook relation.— This relation underestimates the 
exact period for any amplitude, but the difference is almost imperceptible for small angles. 
For larger angles Tq becomes more and more inaccurate for describing the exact period and 
Eq. (1) can be used to obtain a numerical solution. 

Alternatively, an integral expression for the exact pendulum period may be derived from 
energy considerations, without a detailed discussion of differential equations. If we take the 
zero of potential energy at the lowest point of the trajectory (see Fig. 1) and choose for 
simplicity the initial conditions as 6'(0) = +6'o and d9/dt{Q) = 0, we have^ 

1 d6 

mgL{l - cos6'o) = -mL^{—f + mgL{l - cos6'). (2) 
The solution for dO/dt is 



^ = ±y ^(cos6' - cos6'o), (3) 

where the + (— ) sign is for counter-clockwise (clockwise) motion. Now, integrating dO/dt 
from 9q to (thus choosing the — sign in Eq. (3)), which corresponds to a time equal to one 
quarter of the exact period T, we have 

de. (4) 




g Jo \J cos Q — cos ^0 

The definite integral in Eq. cannot be expressed in terms of elementary functions. ^'^ Note 
that the numerical evaluation of the period using Eq. ^ is not straightforward because the 
integrand has a vertical asymptote Q = Q^^ which makes the integral improper.-^ This 
difficulty can be circumvented by substituting cos^^ by 1 — 2 sin^(0/2) and making a change 
of variables given implicitly by sin(y9 = sin(^/2)/ sin(6'o/2). In this way, Eq. (4) becomes 



9 Jo ^/l — k"^ sin^ Lp 



dip, (5) 



where k = sin(^o/2). The definite integral is K{k), the complete elliptic integral of the first 
kind, which is not improper because k < 1 for |^^o| < tt. 

It is not difficult to numerically evaluate T for a given amplitude. The relative error 
made in approximating T by Tq, where T = 2TQ/TTK{k), isf^ 

To-T Tx 

-^ = mk)~ 

Our proposed approximation for the pendulum period is based on the observation that 



k) = ^/l — k'^ sin^ (f is a smooth function of ip, whose concavity changes from downward 
to upward at a point near the middle of the interval of integration, that is, < < 
7t/2. As shown in Fig. 2, this change occurs for any 6*0 between a.nd7T/2M We use the 
points (0,1) and (7r/2,a) for a linear interpolation, where a = f{(p = 71/2, k = sin6'o/2) = 
a/1 — (sin 6q/2Y = cos 6^0/2, and approximate /(</?, k) by 

/(¥.,^o)^l--(l-aV. (7) 

71 

We substitute Eq. ((7j) in the denominator of the integrand in Eq. (0) and find 
Finally, by substituting Eq. (jHJ in Eq. (5), we found 



-lof> 




, = -r.^. (9) 

g 1 — a 1 — a 



Note that Ina < and hence Tiog > for |^^o| < ^- The relative error in the logarithmic 
formula in Eq. (9) is given hy^ 



Tiog — T_ 71 (—Ina) 



T 2K{k) 1-a 



- 1- (10) 



III. COMPARISON WITH OTHER APPROXIMATIONS 

We compare the accuracy of our approximation for the pendulum period in Eq. (jU)) in 
representing the exact period to that of other known approximations for amplitudes less 
than or equal to 71/2}"^ The relative errors found by approximating the exact period by Tq 
and Tiog, as well as by other formulas, are depicted in Fig. 3, where it is easily seen that all 
approximations present the same general behavior: For small amplitudes their corresponding 
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error curves tend to zero and for larger amplitudes the curves go up monotonically, reflecting 
the increase of the relative error with the amplitude obtained with all known approximation 
formulas. However, note that the rate at which the error increases is different for each curve, 
the smaller rate being the desired one for a good approximation. In this sense, the small 
angle approximation T ^ Tq exhibits the worst behavior because its error becomes greater 
than 0.1% (0.5%) above an amplitude as low as 7° (16°). 

The second-order approximation found by Bernoulli in 1749 from a perturbative analysis 
of Eq. (5), perhaps the most famous formula for the large angle period, is^ 

^2 = To(l + ^). (11) 

As seen from the short-dashed line in Fig. 3, it leads to an error that increases rapidly, 
overcoming the 0.1% (0.5%) level for amplitudes above 41° (60°). The addition of more 
terms improves the accuracy of T2.— 

More recently, other approximation formulas have been proposed. Among them, the 
Kidd-Fogg formula has attracted much interest due to its simphcity." It is given by 

Tkf = To^^==- (12) 
a/cos(6'o/2) 

The dash-dotted line in Fig. 3 represents the error for T^f- The error is greater than 0.1% 
only for amplitudes 9q > 57° and reaches 0.8% for 9q = 90°. Thus, it is not accurate enough 
for interpreting the experimental data for very large-angle amplitudes, contrary to the claim 
of Millet.^ 

Another formula for the period arises when an interpolation-like linearization is made 
directly in Eq. (1).— The resulting expression is 

T„.T„«-=". (13) 



which has an error greater than 0.1% only for 9o > 69° (see the thin solid curve in Fig. 3). 
However, the error reaches ~ 0.4% for 9q = 90°, which is four times the typical experimental 
error (0.1%). 

The error using Eq. Q for the period (see the thick solid line in Fig. 3) remains below all 
other error curves for any 9q and is greater than 0.1% only for amplitudes greater than 74°. 
Moreover, it increases slowly, reaching only 0.2% for 9o = 86°. Therefore, our logarithmic 
formula is the better approximation for the exact pendulum period in the sense it yields the 



smaller relative error and also the smaller rate for the relative error increasing in the range 
of amplitudes studied here. 

IV. EXPERIMENT AND RESULTS 

Reliable data for large-angle pendulum periods were obtained by Fulcher and Davis^ using 
a pendulum made with piano wire (measuring two successive swings) and by Curtis'*-^ who 
determined the period as the average of ten successive periods for each initial amplitude. 
Both papers are good examples of accurate period measurements made with an ordinary 
stopwatch. The measurement of the time interval for n successive periods is a good strategy 
for oscillations in the small angle regime, where the amplitude does not change significantly 
from a swing to the next, but not for large-angle oscillations, because the period decreases 
considerably due to air friction. This behavior is confirmed in Fig. 4, where the period T 
in units of Tq is plotted as a function of 9q. In Fig. 4 the curves for each approximation 
formula discussed in Sec. Illll are plotted. Experimental data taken from Refs. \^ andQ and 
the measurements taken by us in a more sophisticated experimenti^ii^ are also shown. The 
experimental data for amplitudes greater than 35° clearly reveal a systematic overestimation 
for the period due to air damping. 

In our experiment both the time-keeping and position detection were done automatically 
to reduce the instrumental error to milliseconds, which is much less than the error in time 
keeping when a common stopwatch is used (of the order of 0.2 s, the average human reaction 
time). We measured the pendulum period by measuring the time interval between two suc- 
cessive passages of the pendulum over the lower point of its circular path, which corresponds 
to T/2. The measurement was based on the variation of the electrical resistance of a light- 
dependent resistor during the passing of the pendulum's bob through the path of the light 
from a laser.^^ Electronic circuitry is needed for converting the analogue signal generated in 
the light-dependent resistor when the pendulum's bob cuts the light's path to TTL compat- 
ible digital voltage, so that the microprocessor can understand the change in current. The 
details of the design/operation of the circuitry and the microprocessor program required for 
measuring the time interval between successive interruptions in the light-dependent resistor 
are in Ref . Q. 

We devoted much attention to the reduction of the air resistance on the motion of the 
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pendulum bob by choosing suitable materials and parameters for the pendulum. We used 
lead as the bob material due to its high density in comparison to other inexpensive metals, 
which gives a small size and large weight for the bob. We found that a cylinder is preferable 
to a sphere because it allows for a better localization of the center of mass, which is needed 
for measuring L accurately. The cylindrical shape also yields a reduction of air resistance 
by reducing the scattering cross-section, that is, by choosing a diameter much smaller than 
the height of the cylinder. These considerations led us to fabricate a body with a mass of 
0.400 kg. For this massive bob we verified that cords made of nylon, a commonly used mate- 
rial, are inadequate because they stretch considerably for large angle oscillations and cause 
undesirable vibrations. The more convenient material taking into account low elasticity. 



lightness (see Ref. |20| for the importance of this factor), price, and availability, seems to be 
cotton. We used a common sewing thread as the pendulum cord. We also investigated what 
cord length would give the best experimental results for large angles.— After comparing 
many lengths for an amplitude of 60°, we choose a length of 1.50 m so that the bob speed 
would be small (because air friction increases with speed, the longer the string length, the 
less the effect of air friction on the period) This length has a period of ~ 2.5 s, which is 
sufficiently small for doing several repetitions of the period measurement for each amplitude 
during a typical one-hour class. These considerations led us to much more accurate experi- 
mental data for the pendulum period for amplitudes less than or equal to 90° as shown in 
Fig. 4. It is seen that our experimental data (black diamonds) are closer to the exact period 
expected in the absence of air resistance (the solid line) than the data in Refs. (crosses) 
and Q (circles). The logarithmic formula in Eq. Q is also in better agreement with the 
experimental data. 
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Figure Captions 




FIG. 1: The pendulum bob is released at rest from a position that forms an angle 6q with the 
vertical and passes at an arbitrary angle 6 (< ^o) with a velocity LdO/dt. Its height depends on 9 
according to L(l — cos 9). 
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FIG. 2: Behavior of the function f{ip, k) = \/l — k"^ sin^ ip for < ip < ir/l rad and for some values 
of ^0 (k = sin(0o/2))- The horizontal and vertical dashed lines are for f(ip,k) = 1 and p = tt/2 
rad, respectively. The dash-dotted lines are the linear interpolation in Eq. for = ^r/G, tt/4, 
7r/3, and 7r/2 rad. 
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FIG. 3: Comparison of the relative errors for the various approximations discussed in the text for 
the period. All curves increase monotonically with ^o- The horizontal dashed line marks the 0.1% 
level. The small angle approximation (T ~ Tq) yields an error that is greater than 0.1% for > 7° 
and reaches 15.3% for > 90°. The thick solid line is for Eq. 0. Note that it remains below all 
other curves for 0° < 6*0 < 90°. 
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FIG. 4: Comparison of the ratio T/Tq for the approximation formulas discussed in the text and 
experimental data. The dotted curve is for the Bernoulli formula, Eq. (11). The dash-dotted curve 
is for the Kidd-Fogg formula, Eq. (12). The dashed line is for our logarithmic formula, Eq. ©. 
The solid line is the curve for the exact period, found numerical integration of K{k). The 
experimental data were taken from Ref. |j (+) and Ref. 
experimental data. 



1^ (o), and the black diamonds are our 
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